Non-Fermi liquid signatures in the Hubbard Model due to van Hove singularities 
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When a van-Hove singularity is located in the vicinity of the Fermi level, the electronic scattering 
rate acquires a non-analytic contribution. This invalidates basic assumptions of Fermi liquid theory 
and within treatments based on perturbation theory leads to a non-Fermi liquid self-energy and 
transport properties. Such anomalies are shown to also occur in the strongly correlated metallic 
state within dynamical mean-field theory. We consider the Hubbard model on a two-dimensional 
square lattice with nearest and next-nearest neighbor hopping within the single-site dynamical 
mean-field theory. At temperatures on the order of the low-energy scale To an unusual maximum 
emerges in the imaginary part of the self-energy which is renormalized towards the Fermi level 
for finite doping. At zero temperature this double-well structure is suppressed, but an anomalous 
energy dependence of the self-energy remains. For the frustrated Hubbard model on the square 
lattice with next-nearest neighbor hopping, the presence of the van Hove singularity changes the 
asymptotic low temperature behavior of the resistivity from a Fermi liquid to non-Fermi liquid 
dependency as function of doping. The results of this work are discussed regarding their relevance 
for high-temperature cuprate superconductors. 
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I. INTRODUCTION 

The question whether interacting fermions form a Lan- 
dau Fermi liquid state at low temperatures or rather 
transform into something commonly referred to as sin- 
gular or no n- Fermi liquid is a very fundamental one. A 
theoretical description requires at least a qualitative un- 
derstanding of the low temperature phase in order to 
address the relevant degrees of freedom in that phase. 
But for strongly interacting electronic systems this in- 
formation is usually lacking, and the low-energy degrees 
of freedom might be very different from free electrons. 
However, the Fermi liquid picture is often appropriate 
at low temperatures and energies. Reasons for this can 
be found in the stability of the Fermi liquid fixed point 
of the renormalization group analysis^ or phase space 
arguments.^ In both approaches, the analyticity of the 
scattering vertices is crucial for the arguments. 

In situations, where the Fermi liquid description fails, 
the origin for this is linked to low-dimensionality of the 
system, long-range interactions, nested Fermi surfaces or 
the proximity to a quantum critical point, just to name 
a few (see, e.g. Refs. 3-6). The quasiparticle concept 
might still be valid, but the scattering rate is strongly 
enhanced and shows an unusual, i.e. non-quadratic, tem- 
perature and frequency dependence. Alternatively the 
quasiparticle concept may also break down completely. 

In this work we will focus on the single-band Hubbard 
model which is a minimal model for strongly correlated 
electron systems, and obtain a non-perturbative approx- 
imative solution by means of the dynamical mean-field 
theory (DMFT)."^^^ Within DMFT, where only local cor- 
relations are included, the generic low temperature state 
of the metallic Hubbard model is a Fermi liquid. *"^^ This 
is expected, since the low temperature phase of the spin- 
^ single impurity Anderson model — the model onto 
which the Hubbard model is mapped within DMFT — 



is a local Fermi liquid under the assumption of a "well 
behaved" medium. ^^ The Fermi liquid reveals itself via 
a low-energy many-body bandstructure forming around 
the Fermi level at low temperatures. ^^'^^ Due to the ne- 
glect of momentum dependent correlations in DMFT the 
heavy quasiparticles possess the non-interacting Fermi 
surface.^ A characteristic low-energy scale Tq is associ- 
ated with this lattice version of the Kondo effect and 
marks the temperature scale for local-moment screening 
and the emergence of coherent quasiparticles. 

On the other hand it is known from perturbation 
theory and renormalization group treatments, that a 
van Hove singularity of the non-interacting density of 
states (DOS) in the vicinity of the Fermi level lead 
to a marginal^^ or non-Fermi liquid form of the self- 
energy. -'^^"^^ This should be contrasted to nested Fermi 
liquids,^ where due to the nesting property of the Fermi 
surface the phase space volume for low-energy scattering 
is also strongly enhanced and unusual low-energy prop- 
erties emerge as well.^'^^ 

Van Hove singularities are a consequence of maxima 
or saddle points in the dispersion relation tk- For peri- 
odic energy bands in a crystal a certain number of these 
van Hove singularities must occur on topological grounds 
(see, e.g. Ref. 33). For example, all of the cubic lattices 
have noticeable van Hove singularities in space dimen- 
sions d ^ 3, typically logarithmic divergences or square- 
root cusps. 

Anomalous low-energy behavior can already be antici- 
pated by recognizing, that the usual arguments of micro- 
scopic Fermi liquid theory break down if a maximum or 
saddle point of the dispersion relation tk is found at the 
Fermi level. Then, tk cannot be approximated by a linear 
dispersion, i.e. tk 76 -^|Afc|, with Afc = k — kp and kp 
the Fermi wave vector. Around these points a quadratic 
(maximum) or hyperbolic (saddle-point) functional de- 
pendence tk ^ Afc^ ± Afcn results, where k±^ and fcy are 



two linear independent directions in momentum space. 
The arguments leading to a quadratic Fermi liquid en- 
ergy and temperature dependence of the scattering rate, 



i.e. ImS'^ {uj - iO+) 



-TT T , are thus not applicable. 



Instead, the phase space volume for scattering is strongly 
enhanced and anomalous T- and w-dependencies may re- 
sult. If the singularity occurs in the vicinity of the Fermi 
level, strong modifications of the low-energy and tem- 
perature properties are still expected due to a pinning 
efrect.26.34-38 

In the present work we show that the non-Fcrmi 
liquid signatures which were encountered in previous 
studies^*"^^ of the self-energy and transport proper- 
ties can also occur in the paramagnetic phase of a 
strongly correlated metal within DMFT. We use the 
two-dimensional square lattice with nearest and next- 
nearest neighbor hopping, t and t' respectively. For 
—0.5 < t' /t < 0.5 the dispersion relation tk has sad- 
dles at the X points in the Brillouin zone {kx = 
{(0, ±7r), (±7r, 0)}). This leads to a logarithmic diver- 
gence in the non-interacting DOS, which is located at 
the Fermi level to = for t' = and moved to finite 
energies for t' ^0. 

The resulting anomalies consist of an unusual maxi- 
mum in the imaginary part of the self-energy IraT,^ {to — 
iO+) in the vicinity of the Fermi level at low but fi- 
nite temperatures. The concept of quasiparticles remains 
valid but they are strongly scattered at the flat parts of 
the dispersion relation and the Fermi liquid formation is 
disturbed. The anomalous double well structure persists 
over a range of doping and next- nearest neighbor hopping 
t'. At zero temperature and for half-filled square lattice 
the double-well structure is removed and the self-energy 
vanishes at the Fermi level. Even though the self-energy 
exhibits a non- Fermi liquid energy dependence, the quasi- 
particle weight is still finite. The pinning of the van Hove 
singularity to the Fermi level, which has previously been 
found, ^^'^^''^^"^^ is explicitly observed in the correlated 
metal. 

Qualitatively similar non-Fermi liquid signatures are 
also observed in more advanced theories which include 
nonlocal correlations. ^^"''^ These features are usually at- 
tributed to the presence of strong nonlocal antiferromag- 
netic correlations. However, in the present approach non- 
local correlations are not included and we argue that an 
interpretation exclusively in terms of nonlocal correla- 
tions is too simplified. 

The paper is organized as follows. In Section II we 
shortly present the Hubbard Hamiltonian and the ba- 
sic equations of DMFT. We elucidate the role of the 
non-analyticities in the non-interacting DOS by means 
of a simplified model DOS for which analytical insights 
are obtained in Sect. HI. Self-consistent DMFT calcu- 
lations for singular and non-singular model DOS reveal 
the non-analyticity to be the cause for the observed non- 
Fermi liquid signatures. In Sect. IV we present finite and 
zero temperature calculations for the spectral function 
and self-energy for a square lattice without next-nearest 



neighbor hopping t' = 0. The anomalies are shown to 
occur at temperatures on the order of the low energy 
scale To and decrease with temperature. At T = a 
non-quadratic energy dependence in the imaginary part 
of the self-energy remains. The case of the square lattice 
with next- nearest neighbor hopping t' = —0.2t is dis- 
cussed in Sect. V. In addition to the spectral function 
and the self-energy, we also study the temperature depen- 
dence of the zero frequency quasiparticle scattering rate 
and the resistivity. We discuss in Sect. VI the relevance 
of the present findings to the cuprate superconductors 
before we conclude the paper. 



II. MODEL AND METHOD 

We study the single band Hubbard model 

^ = Y1 *yCcj<T + Yl e c^Cia + UY^ n,^h,i , (1) 

ij.a i a i 

where the operator Cia (cj^.) annihilates (creates) an elec- 
tron in a localized Wannier orbital at lattice site i with 
spin (T, hju = cj^Cia represents the number operator for 
electrons, e is the local ionic level position at each site, 
and U is the local matrix element of the Coulomb repul- 
sion. 

In the following an external magnetic field is not in- 
cluded and all spin up and down quantities are equal. 
The spin label of Green's functions and self-energies will 
be suppressed. We focus on the paramagnetic phase of 
the model and exclude possible phase transitions to or- 
dered phases such as superconductivity or magnetism. 

All information on the lattice is encoded in the one- 
particle hopping amplitude tij. Its Fourier-transform 
gives the dispersion relation 



tk = — 2<(cos kx + cos ky 



' At' cos kx cos ky (2) 

where we set the lattice constant a — 1 and already speci- 
fied to the two-dimensional (2d) square lattice with near- 
est and next-nearest neighbor hopping, t and t' , respec- 
tively. 

The central quantity for the discussion of Fermi liq- 
uid properties is the correlation self-energy YP {z). The 
assumption that it does not depend on momentum k con- 
stitutes the major approximation of DMFT. The lattice 
Green's function of the Hubbard model is then given by 



G{tk, z) 



z-tk~T.U{z) 



(3) 



which only depends on k via the dispersion relation 
tk- The latter also determines the non- interacting DOS 
Po(w) = jjY.k5{uj-tk). 

Within DMFT the lattice model is mapped onto an 
effective impurity model, where the medium r(z) of the 
impurity model cannot be chosen arbitrarily, but has to 



be determined from the self-consistency conditions 

r(z) 



r(z) 



i + G(z)r(z) 



N^G{z)-^-tk J G(z)-^-x 



G{z) 



1 



S^(z) 



(4) 
(5) 
(6) 



The functional form of the local scattering matrix T[z) is 
completely determined by the non- interacting DOS. The 
correlation self-energy enters only as an unspecified pa- 
rameter [via G{z)\ and is obtained from the solution of 
the effective impurity model for a given medium T{z). 
The self-consistency cycle is closed by identifying the lo- 
cal Green's function of the lattice with the one of the 
impurity model, 



G'(loc)(^)^ 



1 



z-e-YP{z)-T{z) 

1 Y^ Git^, z) = G{z) + G{z) T{z) G{z 



(7) 



TV 



by means of a simplified model DOS of the form 

1 + 



cusp 



{u) 



2aW 



1 



H 
W 



a> 



(8) 



with -W < u <W. For a < 2 this DOS has a non- 
analytic cusp at oj = which turns into a logarithmic 

divergence for a = 0, Pa'^o(w) = jPk ^"^ ( tBt ) ' '^^ile it is 
smooth for a > 2 (see inset of Fig. 1). The case a = 
mimics the logarithmic divergence of the 2d simple cubic 
lattice. 

The explicit form (8) allows for the analytic calculation 
of the effective medium r(z) as shown in the Appendix. 
We state the result for three values of a 



ra=o(2) 



1 



2W 



G{z) U2[-G[z)W] - U2[G{z)W] 



Ta = l{z) 



G{z) 



(9) 



(10) 



W 



G{z)W 



In 



1 - {G{z)WY + 2atanh[G(z)VK] 



For the finite temperature calculations we use the en- 
hanced non-crossing approximation^^' ^^'''^"^^ (ENCA) as 
the impurity solver. There are no adjustable parame- 
ters and the method works directly on the real frequency 
axis. This is crucial for the observation of the small low- 
energy anomalies in the self-energy. However, it cannot 
be employed at too low temperatures T <^To due to the 
NCA pathology ''^ Within the ENCA the pathologic be- 
havior is considerably improved compared to NCA.^° It 
can be reliably used down to temperatures T « Tq/10 
for moderate to large values of U and not too far away 
from half-filling. 

At zero temperature we use the numerical renormal- 
ization group^^"^^ (NRG) instead. For the NRG spectral 
functions we use the complete Fock-space^^'^^ algorithm 
of Ref. 56 with a discretization parameter A = 1.7 and 
keep approximately 1700 states in each NRG-iteration 
step. We average over eight different discretizations^^ of 
the conduction band to minimize the errors. 

We also performed finite temperature DMFT(NRG) 
calculations which were in accord with the 
DMFT(ENCA) results and confirmed the existence 
of the anomalous structures in the self-energy. However, 
the anomalous features occur at energies and tempera- 
tures where the results depend on the actual values of 
NRG broadening parameters. ^^'^* Therefore, we prefer 
the ENCA for the finite temperature calculations. 



III. MODEL DENSITY OF STATES 

Before presenting DMFT results for the Hubbard 
model on cubic lattices we elucidate the basic mechanism 



^W = ^ 



G(z) 



(11) 



{2/3)W 



G(z)W 



1 - 



(G{z)WY 



atanh[G(z)W^] 



where Li2(2;) is the Dilogarithm. The crucial difference 
in the analytic structure of the three results is the ap- 
pearance of non-analytic logarithms for a = and a = 1 
which are absent for a = 2. 

Suppose the system is about to form a Fermi liquid 
state as it is expected within DMFT for the Hubbard 
model at low T . Then, the correlation self-energy has 
the form, 

u-e-YF{LO- iO+) = - (cj - e - i7 [(Trr)^ -H cj^] ) , 

(12) 

where T is the temperature, 7 a measure for the resid- 
ual quasiparticle scattering, and Z the quasiparticle 
weight. The renormalized single-particle energy e = 
e -I- ReS'^(— iO+) determines the Fermi wave vector via 
e = tk and vanishes in symmetric situations. 

Using the form (12) in the expressions Eqs. (9) to (11), 
the imaginary part of F develops a cusp-like minimum at 
the Fermi level for small 7 and a < 2.^^ Figure 1 illus- 
trates this for the three different values of the exponent 
a = 0, 1, and 2. In the limit 7 — >■ the effective media 
for a ~ 1 and q = 2 approach the values ImFa^i = I/tt 
and ImrQ,=2 = 4/(37r), respectively, while for a = the 
effective medium eventually reaches zero (see Eqs. (A. 8) 
to (A. 10)). Thus, for a = 0, we end up with a soft-gap ef- 
fective medium for which Kondo screening and the local 
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FIG. 1: Imaginary parts of the effective media r(a; — iO"'') 
[see Eqs. (9), (10) and (11)] as functions of energy u for the 
cases a = 0,1,2 and two values of the residual interaction 
7. The inset shows the model DOS for the same values of a. 
Vy = Z = 1 and r = ? = was used. 



Fermi liquid behavior can be completely destroyed un- 
der certain conditions. ^°"^^ However, in the present case 
the effective medium vanishes only logarithmically at the 
Fermi level, i.e. Imra=o(w-«0+) ^ -^^ [see Eq. (A. 8)], 
and therefore the low temperature properties arc still ex- 
pected to be characterized by the strong coupling fixed 
point as will be shown in Sect. IV D. 

The topology of the Fermi surface does not enter the 
present arguments at any point, only the non-analyticity 
in the non-interacting DOS does. The latter derives from 
flat parts in the dispersion relation, which in the case of 
the 2d square lattice are saddle points at the X points. It 
is therefore clear, that no nesting property of the Fermi 
surface, which is present in the square lattice at half- 
filling, can be responsible for these effects. 



IV. SQUARE LATTICE 




-0.3-0.2-0.1 0.1 0.2 0.3 -0.2 0.2 

FIG. 2: (a) Local spectral function p(cj) for the symmetric 
Hubbard model on a 2d square lattice for Coulomb repul- 
sion f/ = 4 (e = -2), haff bandwidth W = 2 (t ^ 0.5) 
and various temperatures as functions of frequency within 
DMFT(ENCA). The thin dashed curve represents the non- 
interacting DOS po{i^), and the inset shows the imaginary 
part of the effective medium Imr(a; — iO'*'). (b) Imaginary 
(left) and real (right) part of the self-energy E'^(a; — iO+) for 
the same parameters as in (a) and T as indicated. The inset 
of the left panel shows the imaginary part on a wider energy 
interval. The short arrows in the right graphs point at the 
approximate positions of the kinks. 



A. Half-filling 

In all calculations of this work the hopping matrix ele- 
ment t is used as the energy scale and energies are mea- 
sured in units 2t = 1, setting the effective half bandwidth 
W = 2dt = 2. In this section the next-nearest neighbor 
hopping is set to zero, t' = 0. The non-interacting DOS 
of this 2d square lattice has the same logarithmic singu- 
larity at oj = as the model DOS of Eq. (8) for a = 0. 

Results of a DMFT(ENCA) calculation for the half- 
filled n = 1 Hubbard model on a square lattice are shown 
in Fig. 2 for various temperatures. In part (a) the non- 
interacting DOS as well as the fully interacting local spec- 
tral function p{(jj) = —IinG'-^°'^\u; — iO~^) are displayed for 
U = 2W — 4. The well separated Hubbard bands around 
the ionic level positions w ss ztU/2 = ±2 are clearly visi- 
ble. Upon lowering the temperature the lattice version of 



the Kondo-effect leads to the build-up of a pronounced 
many-body resonance at the Fermi level. 

The van Hove singularity in po leads to the unusual 
minimum in the imaginary part of the effective medium 
114^(0; — iQ~^) [see inset of Fig 2(a)] as it was already 
found in the previous section. In contrast to Fig. 1 the 
minimum in Imr(a; — iO"*") is not cusp-like at w = but 
sn4eared out due to the imaginary part of the self-energy 
which is shown in part (b) of the figure. However, the 
non-analytic structure in the lattice scattering matrix, 
Eqs. (5) and (9), is sufficient to produce this minimum 
at low temperatures. 

The self-consistency of DMFT leads to an accompa- 
nying maximum in Iml]^(a; — iO"*") at a; = 0, which is 
ii4 striking contrast to Fermi-liquid behavior ImY,^ {u! — 
iO+) ^ u;"^ + n'^T'^ for small uj and T. The position of 
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FIG. 3: Position of the minimum ojmin in ImE in units of the 
low-energy scale To as a function of the rescaled temperature 
T /To for three different Coulomb repulsions. The inset shows 
the screened local magnetic moment normalized to its high 
temperature value /x^ as a function of temperature for the 
same values of the Coulomb repulsions. The horizontal line 
indicates the value where /^e// reaches 70% of ^^o and the 
vertical dashed lines indicate the low-energy scales: To{U — 
3.5) f» 0.082, To{U = 4) ^ 0.051, and To{U = 4.5) f» 0.023. 



energies as well as the minima at zbwmin are symmetric 
around cj = for half- filling) . For all U, Wmin/^o exhibits 
a dome-shaped curve with both, maximum position and 
height determined by Tq. The physical origin can there- 
fore be directly linked to the emergence of the low-energy 
quasiparticles. At high temperatures the system is dom- 
inated by the incoherent charge and spin excitations and 
the van Hove singularity is concealed. At temperatures 
of the order of Tq the lattice version of the Kondo ef- 
fect leads to the screening of local magnetic moments. 
The break-up of screened magnetic moments gives rise 
to long-lived low-energy quasiparticle excitations which 
manifest themselves in the many-body resonance at the 
Fermi level as observed in the spectral function of e.g. 
Fig. 2(a). But their evolution into the coherent Fermi 
liquid quasiparticle is disturbed by the enhanced scatter- 
ing due to the van Hove singularity. This leads to the 
maximum in ImS^ at the Fermi level. Further lowering 
the temperature the phase space volume for quasiparti- 
cle scattering shrinks and so does the extension of the 



maximum. 



B. Finite doping 



the two minima emerging at finite frequencies next to 
the central maximum is moved toward the Fermi level 
with decreasing temperature so that the anomalous re- 
gion shrinks. 

As a side remark we add here, that we observe addi- 
tional kinks in the real part of the self-energy as shown in 
Fig. 2(b), which are associated with this non- Fermi liquid 
behavior. At energies uj w ±0.2 the known kinks aris- 
ing from the coupling of the quasiparticles to local spin 
fluctuations are observable. ^^'^^ The additional kinks are 
clearly visible at lower energies. In contrast to the kinks 
at higher energy, their position is temperature dependent 
and moves toward the Fermi level with decreasing T. 

Anomalous maxima in ImS'^ are also encountered 
within the functional renormalization group approach to 
the Hubbard model at finite temperatures. '^"''^"'^ Crucial 
for these anomalies to occur is the renormalization of 
the two-particle interaction vertex. Within DMFT all 
local n-particle vertices and their renormalizations are 
included. 

We extract the low-energy scale Tq as the temperature 
where the effective screened local moment as calculated 
from the local magnetic susceptibility via^° /i^jj = ^x"oc^ 
is reduced to 70% of its high temperature value /i^ (see 
inset of Fig. 3). This yields the equivalent result to the 
width of the many-body resonance at the Fermi level in 
the momentum resolved spectral function ImG'(ife , lu — 
iO~^)/Tr (not shown). 

The temperature evolution of the anomalous double- 
well structure is governed by the low-energy scale as re- 
vealed by Fig. 3. There, the position of the anomalous 
minimum Wmin in ImE^ as a function of temperature for 
three values of the Coulomb repulsion is shown (the self- 



The van Hove singularity is shifted away from the 
Fermi energy uj = when doping the system, n = 1 — S. 
The corresponding anomalous minimum in the effective 
medium ImF and the maximum in ImE^ are too moved 
to finite energies. This can be seen in Fig. 4, where 
the spectral function, the self-energy, and the effective 
medium are shown for the Hubbard model on a 2d square 
lattice {f = 0) with t/ = 6 and 5 = 0.02 (fiUing n = 0.98). 

The quasiparticle peak around w = in the spectral 
function, as displayed in the inset of Fig. 4(a), acquires 
a pronounced asymmetry due to the double-well struc- 
ture in ImS'^. The shoulder in the curve for the lowest 
temperature might even be interpreted as a precursor of 
a pseudogap. 

The flat parts in the dispersion relation arc located en- 
ergies away from the Fermi level and the dispersion rela- 
tion can again be approximated linearly in a vary narrow 
region around the Fermi surface. In this region the usual 
arguments of microscopic Fermi liquid theory hold and 
the system forms a regular Fermi liquid at very low tem- 
peratures. This is observed in the self-energy shown in 
Fig. 4(b) where a quadratic minimum forms at w = at 
T ~ 0.007. However, the Fermi liquid parameters will be 
strongly rcnormalized as compared to a situation without 
the van Hove singularity. 



C. Self-consistent treatment with a model DOS 

No true divergence in the non-interacting DOS is nec- 
essary for the anomalies described above to occur. A 
non-analytic cusp as in the model DOS of Eq. (8) for 
< a < 2 is sufficient. This is illustrated in Fig. 5 where 




FIG. 4: (a) Local DOS for the asymmetric Hubbard model 
on a 2d square lattice with U — 6 and filling n — 0.98 for 
various temperatures obtained with the DMFT(ENCA). The 
inset shows the low-energy region, (b) Imaginary part of the 
self-energy around the Fermi level for the same parameters 
as in (a) and temperatures as indicated. The inset shows 
the imaginary part of the effective medium in the low-energy 
region. 



fully self-consistent DMFT calculations for several model 
DOS are considered. The first two plots (a) and (b) are 
for p™*'P(a;) with a = 1 at and away from half-filling, re- 
spectively. Figures (c) and (d) are for asymmetric DOS 
which are constant below the Fermi level and have linear 

(c) and a square-root (d) cusp at w = 0. The imaginary 
parts of the effective media and self-energies show the 
characteristic minima and maxima induced by the non- 
analyticities in pg. For the half-filled cases (a), (c) and 

(d) the anomalies are situated at the Fermi level, while 
for finite hole doping as in (b) these are moved to positive 
energies. 



D. Zero Temperature 

Figure 6(a) displays the spectral function for the Hub- 
bard model on a square lattice with [/ = 4 at zero tem- 
perature and for different doping 6. These spectra were 
obtained with the NRG^^ as impurity solver. 
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FIG. 5: The imaginary parts of the effective medium (up- 
per graphs) and self-energy (lower graphs) of the Hubbard 
model within DMFT(ENCA) for different non-interacting 
DOS, which are shown in the insets. The vertical dashed 
line indicates the position of the Fermi level. The parameters 
are indicated in the plots. 



As the hole-doping is increased, the lower Hubbard 
band moves toward the Fermi level and eventually merges 
with the quasiparticle peak. For very large doping the 
lower Hubbard band is moved above the Fermi level and 
p{u!) resembles the non-interacting DOS as the system 
becomes effectively non-interacting. 

The position of the quasiparticle peak is attracted 
to the Fermi level up to considerable values of doping 
(roughly 6 < 0.2). This pinning becomes more visible in 
the inset, where the position Wmax of the maximum in 
p{Ld) is plotted as a function of doping S ioi U — 4 and 
U — 0. In the interacting case Wmax depends quadrati- 
cally on the doping, Wmax ^ S'^, and is considerably re- 
duced compared to U ~ 0. The many-body correlations 
renormalize the van Hove singularity toward the Fermi 
level especially at small doping i5 < 0.2. 

The imaginary part of the self-energy shown in 
Fig. 6(b) does not display any anomalous double- well 
structure. The logarithmic van Hove singularity in po 
is a rather weak divergence and at zero temperature the 
usual arguments leading to a vanishing self-energy at the 
Fermi level apply. ^"^ However, ImE^ shows an anoma- 
lous frequency dependence at low energies, as it grows 
like ImE'-'(a; — lO"*") ^ jwp/^ which was also previously 
found in perturbation theory.^*^'^^^*^ The upper inset re- 
veals this low energy behavior more clearly in a double 
logarithmic plot. 

The quasiparticle peak at the Fermi level for half-filling 
i5 = shows a logarithmic divergence as was already 
observed recently. ^^ But in contrast to what was stated 
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FIG. 6: (a) DMFT(NRG) spectral function at T = for the 
Hubbard model on a 2d square lattice for f/ = 4 and varying 
hole doping. The inset shows the position of the maximum as 
a function of doping for [7 = 4 and C/ = 0. The dashed line 
is a fit Wmax = 2.25^ . (b) Imaginary part of the self-energy 
at half-filling J = for three values of the Coulomb repulsion 
U. The curves for U = 4.5 and U = 5 were rescaled in order 
to lie on top of the (7 = 4 curves at low to. The dashed line 
is a fit with a function a\uj\ ' . The upper inset shows the 
low energy behavior in a double- logarithmic plot. The lower 
inset shows the quasiparticle weight for [/ = 4 as a function 
of doping. 



in Reference 71, the quasiparticle weight 
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(13) 



0=0 



does not vanish at half-filling in our calculation, as can be 
seen in the inset of Fig. 6(b). (The quasiparticle weight 
in Fig. 3 of Ref. 71 also seems to extrapolate to a finite 
value at half-filling but a contradicting statement is made 
in the text of that reference.) The slope of the real part 
of the self-energy does consequently not diverge at the 
Fermi level. This is in accord with the low frequency 
dependency of the imaginary part ~ Iwp^ from which 
the real part is obtained via Kramers-Kronig relation. 
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FIG. 7: DMFT(ENCA) spectral functions p{io) (left panels) 
and imaginary parts of the self-energy TF (uj — iO^) (right 
panels) for different filling n for the Hubbard model on the 
square lattice with additional next-nearest neighbor hopping 
t' = -0.2i = -0.1 and U = 4.25. The insets in the left panels 
show the imaginary part of the effective medium while the 
insets on the right show a close-up of the low-energy region 
of ImE'^ . The non-interacting DOS is shown in the inset for 
the half- filled case (b). 



ReH^ {oj — iO+) ~ uj[y/uJ{a.cothy/uJ + arccot-^o;) — 2]. 

The system is therefore well characterized by a gener- 
alized Fermi liquid^^"®'^ where ImS'^ shows some anoma- 
lous frequency dependence, but still vanishes at the Fermi 
level and the quasiparticle weight remains finite. 



V. NEXT-NEAREST NEIGHBOR HOPPING 

A. Spectral function and self-energy 

In this section we extend the study to an additional 
next-nearest neighbor hopping t' = — 0.2<. The local 
spectral functions and self-energies displayed in Fig. 7 
for U = 4.25, various fillings and temperatures show pro- 



nounccd particle-hole asymmetries. 

At low temperatures the anomalous minima in ImF 
and corresponding maxima in the self-energy ImS'^ can 
be observed slightly below the Fermi level for filling n = 
1.04, n = 1, and n = 0.98. This correlates with the 
logarithmic van Hove singularity in po shifted below the 
Fermi level, as displayed in the right inset of Fig. 7(b). 
For decreasing filling (increasing hole doping) the van 
Hove singularity is moved toward the Fermi level and 
the anomalous structures in ImE and ImF narrow. 

For the largest hole doping n = 0.92 the imaginary part 
of the self-energy exhibits only one narrow minimum at 
very low temperature. Within DMFT the Fermi surface 
does not change upon increasing the interaction U,^ due 
to the momentum independence of the self-energy. The 
filling at which the van Hove singularity moves across the 
Fermi level can be calculated with the non-interacting 
DOS to be n sa 0.82, which is much less than n = 0.92 
Therefore a double-well structure would be expected but 
is not observed. The reason is found in the lower Hub- 
bard band reaching the Fermi level. A growing number of 
empty lattice sites arc created and disturb the quasiparti- 
cle formation. Additional c^uasiparticle scattering results 
and conceals the effect of the van Hove singularity. 

Due to the pinning of the singularity to the neighbor- 
hood of the Fermi level its influence is strong for all values 
of filling. 

The anomalies correlate with the low-energy scale and 
occur at temperatures following the trend given by To as 
function of filling n. Tq reaches a minimum for n sa 0.98 
to 1 where correlations are strongest and the system is 
closest to the Mott insulator. 



B. Transport 

The anomalous behavior found in the self-energy af- 
fects the temperature dependence of transport proper- 
ties. The static quasiparticle scattering rate at the Fermi 
level Iml]'^(— iO+) is shown in Fig. 8(a) as a function of 
temperature for different filling. 

For a 2d Fermi liquid the expected temperature depen- 
dence of the scattering rate is given by ImE'^(— iO+, T) ^ 
T^ ln(l/T). This behavior is found for low temperatures 
at half-filling [lowest curve in Fig. 8(a)] and above, i.e. 
n > 1 (not shown). Hole-doping induces qualitative 
changes in the low temperature behavior since the van 
Hove singularity is moved toward the Fermi level. For fill- 
ing n = 0.98 and n = 0.96 power-law behavior with an ex- 
ponent close to one is observed, Iml]'^(— zO'*', T) ^ T^-^ . 
For large hole doping, n ~ 0.92, the scattering rate is 
again found to be in accord with usual Fermi liquid vari- 
ation, ~ T^ln(l/T). This is expected to hold for even 
larger hole-doping n < 0.92, but could not be checked 
within the present approach due to restriction of the 
ENCA to relatively high temperatures at large doping. 

We want to emphasis that the van Hove singularity in- 
duces qualitative changes in the asymptotic low temper- 



10^ 



lO'' 



I ml 



W 
B 



10" 
10" 





Xl01)(] 

^..i,.i-»«"-«-"-"'*x*i', 

-■«■-""' 
„.«»^ 


•■^ 0.98 ■«- 


y^ 


0.96 -.■ - 


(a) -~T2hi(a/T) 


0.92 -n- 



0.01 



0.1 



105 


- n = 1 — 






0.98 ".-■ 


_ . ^,„„-i ■ 


10" 


- 0.96 ■■»■■ ^^ 


^ 




0.92 ..».. ~nr^ 


,.»■•••*'"* 


10-' 


.*"* 






^' 


^^ 


10^ 


'^J^-^ 


y^ 300 


10' 


'~T^ j/" 


200 




-T?-*^ 


100 


10" 


" (b) 










0.5 



0.01 



0.1 



FIG. 8; (a) Static quasiparticle scattering rate ImE (— iO"'") 
for the Hubbard with U = 4.25 on a square lattice with next- 
nearest neighbor hopping t' = — 0.2f as function of tempera- 
ture in log-log plot. The curves for different filling are mul- 
tiplied by the factors indicated on the right, (b) Resistiv- 
ity Q (normalized to ^o) for the same parameters as in (a). 
The inset shows q without the scaling factors and in a non- 
logarithmic plot. The dashed lines represent fits with the 
indicated functions in both plots. 



ature behavior of the scattering rate when moved closer 
to the Fermi surface by doping. ImE^(— lO+jT) de- 
creases much slower with temperature for fillings around 
n K, 0.96 when the double-well structure is found in the 
vicinity of the Fermi level. 

With the self-energy at hand one can also calculate 
transport coefficients. No vertex corrections occur in 
DMFT due to the momentum independence of the self- 
energy and two-particle vertex. ^^"^^ The current-current 
correlation function is completely determined by the 
particle-hole bubble. Transport quantities can be ex- 
pressed in terms of a generalized transport lifetime (see, 
e.g., Ref. 8 or 12) 
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where /(w) — I/(e"' -f- I) is the Fermi function. The 



generalized transport lifetime is given by 



-^^H-]^E(|r) [l-G(%,a.-.0+r (15) 

k ^ ^ 



-~^Jdepo{e)[lmG{e,uj-iQ+)y 



(16) 



The function 



«(.)4e(S) *<-'i) ™ 



can be calculated for the cubic lattices and turns out to 
be a smooth function where the singularities of po &i'6 
removed due to the derivative. 

Only in the Fermi liquid regime is Tj;a;(a;) ~ 
1/lmT,^ {u! — iO"*") and the linearized Boltzmann theory 
is recovered 7^ 

The resistivity 
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go 
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(18) 



with Qq = 2Il§ (a is the lattice spacing, e the electronic 
charge) is shown in Fig. 8(b) for different filling. Typical 
characteristics of strongly correlated systems can be ob- 
served for not too large doping. The resistivity g changes 
from insulating (or semi-conducting) behavior at high 
temperatures to metallic behavior at low temperatures 
due to the emergence of (coherent) low-energy quasipar- 
ticles (see inset). ^^ The high-temperature resistivity is 
much larger than the value expected from the loffe-Regel 
condition valid for usual metals, where g is bound by 
a minimal scattering length of the order of the lattice 
spacing. ^^ At intermediate temperatures 0.1 ^ T < 0.5 
linear regimes with increasing slope for decreasing doping 
can be recognized, in accord with earlier studies. *'"'^° 

The log-log plot (main panel) reveals that the low tem- 
perature Fermi liquid behavior g{T) ^ T^ is found for 
n = 1 and above (not shown). '^^ A T^-dependency is also 
observed for n — 0.92. But for the intermediate values 
n = 0.98 and n = 0.96 a qualitative different behavior is 
found. g{T) more closely resembles power-laws with ex- 
ponents 1 and 1.5, respectively. The van Hove singularity 
is close to the Fermi level and the enhanced scattering in- 
creases the resistivity. 

In summary, the interplay between the reduction of 
phase space volume for scattering and the qualitative 
changes in the frequency dependence of the self-energy 
within the temperature window causes the qualitative 
changes in the transport properties. 



VI. RELEVANCE FOR THE THEORY OF 
CUPRATE SUPERCONDUCTORS 

The t-i'-Hubbard model might be suitable to describe 
the low-energy physics of cuprate superconductors. Even 
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FIG. 9: Generalized quasiparticle weight for the Hubbard 
model on a square lattice with next-nearest neighbor hop- 
ping t' = — 0.2f and U = 4.25 as function of temperature for 
three different filling obtained from DMFT(ENCA) calcula- 
tions and Eq. (20). 



though DMFT is non-perturbative and includes nontriv- 
ial local many-body correlations, it represents a rather 
poor approximation for the two-dimensional cuprate lay- 
ers due to the neglect of nonlocal correlations. Therefore, 
we want to point out that we neither try to model the 
cuprate superconductors nor can our results be directly 
transferred to those systems. 

However, even in more appropriate theories like 
clustcr-DMFT''^'' the van Hove singularities associated 
with the lattice structure are retained and therefore the 
findings of this work have some bearings on the pseu- 
dogap and strange metal phase of cuprate superconduc- 
tors. In the following we will elaborate on some of these 
aspects. 

The van Hove singularity provides an additional mech- 
anism for kinks®^'^^'^^"^^ in the real part of the self- 
energy [see Fig. 2(b)]. The location of the kinks is tem- 
perature dependent which should make it possible to dis- 
criminate these type of kinks in experiinent or calcula- 
tions. 

The non-Fermi liquid self-energy induces an unusual 
temperature dependence of the generalized quasiparticle 
weight*^ 



Z{T) = 



1- 



ImS'^'(iLJo) 



1 I r d^ ImSC/(^-iO+) 



(19) 
(20) 



where ojo = t^T is the smallest fermionic Matsubara fre- 
quency (The form (19) is useful when the self-energy is 
obtained only at the imaginary Matsubara frequencies 
like in quantum Monte-Carlo calculations). The gener- 
alized quasiparticle weight as calculated with Eq. (20) is 
shown in Fig. 9 for the Hubbard model on the square 
lattice with next-nearest neighbor hopping as function of 
temperature and three different filling. All curves show a 
change in curvature at the temperature where the anoma- 
lous double- well structure emerges. This behavior closely 
resembles that found in Ref. 83 for a crossover from a 
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marginal Fermi liquid at high temperatures to a Fermi 
liquid at low T. Since the self-energies do not show any 
signs of marginal Fermi liquid behavior (see the previous 
sections) this reveals a difficulty when interpreting imag- 
inary Matsubara data at finite temperatures. Due to the 
integral in Eq. (20) the temperature dependence of Z{T) 
is not determined by the asymptotic low-frequency form 
oilmSP {uj—iQ'^) but by its average spectral weight in the 
temperature window |a;| < ttT. Therefore, the discrimi- 
nation between different characteristic low-energy forms 
of YF {lo — jO"'") given the self-energy at imaginary Mat- 
subara frequencies seems rather delicate at finite temper- 
ature. 

The partial destruction of the Fermi surface (Fermi 
arcs) and the emergence of a pseudogap^^"*^ is of cur- 
rent interest. These phenomena are usually explained to 
arises from the coupling to strong nonlocal antiferromag- 
netic correlations.^ In the present study no nonlocal 
fluctuations are included but instead the van Hove sin- 
gularity produces similar features such as a non-Fermi 
liquid maximum in IthYP {uj — zO"*") and the concomitant 
reduction of spectral weight [see Fig. 4(a)], However, a 
true pseudogap is not observed in our approach indicat- 
ing the importance of non-local fluctuations. 

The existence of Fermi arcs and the asymmetry for mo- 
mentum vectors along the nodal k = A(7r, tt) (0 < A < 1) 
and anti-nodal k = A(0,7r) direction is natural in the 
present scenario. The flat parts in the dispersion relation 
occur at the X point in the Brillouin zone where the low- 
temperature quasiparticles are most strongly affected and 
consequently the Fermi liquid description breaks down 
there first. '^'''^-'^ Also the stability of the induced non- 
Fcrmi liquid behavior up to a doping of the order of 
20% which was encountered in other studies^^ can be 
explained by the pinning of the van Hove singularity to 
the Fermi level. 

Unarguably, nonlocal correlations are vital for the un- 
derstanding of cuprate superconductors. But as the 
above aspects suggest, some of the non-Fermi liquid sig- 
natures can be produced - or at least aided - by the 
presence of a van Hove singularity in the vicinity of the 
Fermi surface. This is especially important for lattices 
without nesting or with frustration, where nonlocal an- 
tiferromagnetic spin-fluctuations are usually suppressed 
and thus might be too weak to produce non-Fermi liquid 
signatures. Therefore, we argue that the interpretation 
solely in terms of collective modes which couple to the 
electronic degrees of freedom is too simplified and the in- 
fluence of the saddle points in the Brillouin zone has to 
be reconsidered. 



VII. CONCLUSIONS 

In the present work we have focused on the influence of 
a van Hove singularity in the non-interaction density of 
states on the low temperature properties of the Hubbard 
model. In order to include non-perturbative correlation 



effects we used the DMFT to calculate the spectral func- 
tion and self-energy. As impurity solvers we employed the 
enhanced non-crossing approximation for finite tempera- 
tures T > and the numerical renormalization group for 
T = 0. Both methods yield dynamic quantities directly 
on the real frequency axis, which avoids the inaccura- 
cies connected with a numerical analytic continuation of 
imaginary time data. 

The van Hove singularity causes profound changes of 
the low-energy Fermi liquid properties usually encoun- 
tered within DMFT for the Hubbard model. For the 
strongly correlated metal close to the Mott Hubbard 
metal-insulator transition the imaginary part of the self- 
energy develops an unusual double- well structure at finite 
temperatures. This anomalous structure appears at tem- 
peratures on the order of the many-body scale Tq which 
also determines its energy spread. It originates from an 
enhanced scattering of the low-energy quasiparticles at 
the saddle points in the dispersion relation associated 
with the van Hove singularity. 

Using a model DOS, we have shown analytically that 
the non-Fermi liquid signature in the quasiparticle self- 
energy is directly linked to non-analytic logarithmic con- 
tributions to the lattice scattering matrix. As a conse- 
quence the medium for the effective impurity develops a 
dip at the van Hove singularity. In case of a logarithmi- 
cally diverging non-interacting DOS this dip produces a 
soft-gap in the effective medium at T = 0. 

At zero temperature the interacting spectral function 
for a square lattice exhibits a logarithmic divergence for 
half-filling. The system is well characterized by a general- 
ized Fermi liquid where the quasiparticle weight remains 
finite and the imaginary part of the self-energy vanishes 
at the Fermi level. But the frequency dependent scat- 
tering rate increases with a non-Fermi liquid exponent 
« 1.5 ^ 2. Upon doping a pinning of the van Hove sin- 
gularity to the Fermi level is observed. Its signature in 
the spectral function is always renormalized toward the 
Fermi level and stays close to it up to rather large doping 
(5<0.2. 

For the square lattice with a finite next-nearest neigh- 
bor hopping t' = — 0.2i the van Hove singularity in the 
non-interacting DOS is located below the Fermi level and 
moved toward it upon increasing hole-doping. This trend 
is also observed in the anomalous maximum in the self- 
energy. 

In case of half-filling and electron doping, as well as for 
large hole doping, the zero-frequency quasiparticle scat- 
tering rate exhibits a low temperature dependence which 
is consistent with that of a two-dimensional Fermi liquid. 
In between these fillings the decrease is much slower and 
follows a power law with exponent close to one. 

A similar signature of the van Hove singularity was 
found in the resistivity. A Fermi liquid T^ dependence 
emerged for n > 1 and n < 0.92, while for 0.92 ^ n < 1 
the resistivity was enhanced and a decrease with expo- 
nents less than two was observed. 

The findings of this work bear some implications on 
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theories for cuprate superconductors. Even thought no 
nonlocal fluctuations were included, we still produced 
qualitative features usually only obtained within more 
advanced theories where those correlations are incorpo- 
rated. This raises the questions concerning the origin and 
physical mechanism behind such features and we think 
the role of the van Hove singularity in connection with 
strong correlations should be further explored in the fu- 
ture. 
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Appendix: Details of the model DOS calculation 

The scattering matrix can be calculated analytically 
via Eq. (5) for any G (omitting the arguments). 



T = 



-w G^^ — X 



a>{) 1 +a 1 



a G 
{GWf 
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3-Ha 

0=0 1 
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3 -I- a 5 -I- a 
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{GWf 
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2 ' 2 
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2GW 



(A.l) 
(A.2) 



(A.3) 



where 2F1 is the Gauss's hypergeometric function and Li2 
the Dilogarithm. Equation (4) then yields the effective 
medium 



J 



G (1 + Q.)atanh(GVK) - GW 

a=0 l_ ^ 2W_ 

G Li2(-GVK)-Li2(Gl^) 



Wa 



'^{Gwf 2F1 (1, 3±^; 5±£i; i^GWy) 



(A.4) 



(A.5) 



For a = 0, Eq. (9) already constitutes the final expres- 
sion, while for the a = 1, 2 equation (A.4) can be further 
simplified. 

For a = 1 the hypergeometric function reduces to 



2Fi(l,2;3;z^)=-^[z^-Hln(l-z^)] , (A.6) 

leading to the effective medium displayed in Eq. (10). In 
case of a = 2 

^Fi(l,^;^;zA ^ -~[3z + z' - 3atanh{z)] , 

(A.7) 

yields the effective medium of Eq. (11). 

Assuming a Fermi liquid as in Eq. (12) and taking the 
limit 7 — > the effective media take the approximate 



low-energy forms, a;<VF=l (Z = 1,T = 0, e = 0) 



imr..o(c^-.o^)^-^ . .;"'"' ,, 

7r2^ + (ln|w|)2 



(A., 



ImrQ,=i(a; — iO^) 



ImrQ=2('^ ~ *0^ 



7r(l-|c.|) 



7r2(l-|a;|)2+4tj2(l_ln|tj|)2 

(A.9) 

" ^-'^' . (A.IO) 



3^ + (4-f)o 



For a = the effective medium reaches zero at zero en- 
ergy. For q; = 1 it saturates at a value I/tt but the 
non-analytic cusp still remains, while for a = 2 the cusp 
is removed and ImrQ=2(— *0+) = 4/(37r). 
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